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I .  INTRODUCTION 


The  aerospace  developments  of  the  last  decade  stimulated 
a  search  for  means  and  methods  of  determining  relevant  flow 
field  parameters.  These  were  important  for  the  understanding 
of  the  phenomena  encountered  in  space  travel  and  reentry 
as  pertaining  not  only  to  the  fluid  dynamic  phase  of  the 
problem,  but  also  to  problems  of  communication  and  discrimination. 
Among  the  important  parameters  to  be  determined  relevant  to 
the  problem  at  hand  are  the  flow  field  configuration,  neutral 
and  ion  densities,  temperature,  etc.  To  gain  an  understanding 
of  phenomena  observed,  it  was  necessary , simultaneously  with 
theoretical  modeling, to  obtain  experimental  data, not  only  from 
on-board  measurements,  but  also  from  ground  simulation 
facilities.  To  obtain  these  data, new  diagnostic  instrumentation 
and  techniques  had  to  be  developed.  In  the  development  of 
these  instruments  one  was  aiming  at  developing  an  instrument 
which  should  combine  simplicity,  economy  and  ease  of  operation. 
Since  in  most  cases  under  consideration  the  flow  was  slightly 
ionized,  a  natural  choice  for  a  diagnostic  instrument  was 
the  electrostatic  probe. 

As  is  well-known,  the  electrostatic  probe  has  been 
treated  in  innumerable  publications  over  the  last  half  a 
century  (Refs.  1,2,3).  Although  in  principle  the  electro¬ 
static  probe  is  the  most  simple  instrument  available,  the 
interpretation  of  the  measured  collected  currents  is  very 
complex.  Depending  on  the  regime  of  operation  of  the  probe 
and  its  geometry,  different  theoretical  as  well  as  empirical 
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techniques  for  the  interpretation  of  the  collected  currents 
in  terms  of  the  corresponding  ionized  particle  densities 
emerged  (Refs.  4,5,6).  These  techniques  were  applied  in 
conjunction  with  other  techniques  and  found  to  be  very 
satisfactory  when  applied  properly. 

The  application  of  the  electrostatic  probe  as  a  field 
tracing  diagnostic  tool  has  been  proposed.  This  diagnostic 
technique  appeared  to  have  the  necessary  merits  for  this 
purpose  and  particularly  in  highly  rarefied  flow  regimes 
where  other  techniques  and  particularly  optical  techniques 
became  inapplicable.  In  utilizing  this  field  tracing 
technique  in  ground  based  simulation  facilities  (shock 
tunnels,  ballistic  ranges),  two  conditions  have  to  be  satisfied; 
namely,  one  must  know  the  flow  regime  (Xn_n)  beforehand  and 
the  flow  must  have  a  certain  degree  of  ionization.  As 
indicated  in  some  of  the  references  (Ref.  7  and  8),  flow  field 
configurations  were  obtained  using  the  electrostatic  probe 
as  the  diagnostic  tool.  The  technique,  however,  as  pointed 
out  in  Ref.  7  can  only  be  applied  in  nondiffusive  regions  if 
the  ionized  particle  densities  are  to  be  interpreted  in  terms 
of  the  neutral  density  of  the  flow  field.  Since  the  flow  fields  in  the 
above  references  consisted  of  diffusive  and  nondiffusive 
regions,  and  in  view  of  the  restriction  placed  on  this 
technique,  it  became  clear  that  these  field  configurations 
could  be  interpreted  on  a  qualitative  basis  only.  In  order 
to  obtain  a  better  understanding  of  the  phenomena  occurring 
in  these  flow  fields,  it  became  necessary  to  determine  flow  field 
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configurations  utilizing  a  different  diagnostic  technique. 

One  of  the  most  attractive  tools  for  this  purpose  appeared 
to  be  the  electron  beam  density  probe.  As  discussed  in  the 
literature,  this  technique  has  been  utilized  extensively  for 
this  purpose  (Refs.  9-11).  It  was,  therefore,  felt  that  this 
diagnostic  tool  in  conjunction  with  the  electrostatic  probe 
tracing  technique  could  be  utilized  in  a  complimentary 
fashion.  It  should  be  pointed  out  that  the  electron  beam 
density  probe  in  conjunction  with  the  proper  calibration 
technique  is  capable  of  providing  absolute  density  measurements. 
It  would,  therefore,  be  possible,  utilizing  the  electron  beam, 
to  define  the  region  of  applicability  of  the  electrostatic 
probe  tracing  technique.  In  view  of  the  above,  an  experimental 
program  was  initiated.  It  should  be  mentioned  here  that 
consistent  with  range  of  applicability  of  these  techniques 
the  experimental  program  was  designed  for  the  slightly 
ionized  rarefied  hypersonic  flow  of  the  PIB  hypersonic  shock 
tunnel.  This  program  consisted  of  obtaining  flow  field 
configurations  of  models  where  diffusive  as  well  as  nondiffusive 
regions  existed.  In  addition,  flow  fields  of  a  highly 
complex  nature  were  traced.  The  investigated  flow  fields 
consisted  of  the  near  wake  of  a  pointed  10°  half  angle 
cone  at  zero  angle  of  attack,  and  the  disturbed  region  about 
a  sharp  leading  edge  flat  plate.  Both  ion  and  neutral  density 
profiles  were  obtained.  Comparison  between  the  flow  fields 
as  obtained  utilizing  these  two  diagnostic  techniques  were 
made  for  both  models.  In  addition,  the  obtained  measurements 
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over-  the  plate  were  compared  with  calculated  values  using 
the  Rudman-Rubin  solution  for  the  hypersonic  flow  about  a 
sharp  leading  edge  flat  plate.  This  solution  was  modified 
to  include  ionized  particle  behavior  in  the  flow. 

The  authors  would  like  to  express  their  thanks  to 
Dr.  M.H.  Bloom,  Dean  of  Engineering  and  Director  of  Gas 
Dynamics  Research,  for  his  encouragement,  stimulating  dis¬ 
cussions  and  support  in  the  course  of  this  work. 
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II.  THEORY  OF  SLIGHTLY  IONIZED  HYPERSONIC  W 
ABOUT  A  SHARP  LEADING  EDGE  FLAT  PLATE 

The  slightly  ionized  hypersonic  flow  about  a  sharp 
leading  edge  flat  plate  is  discussed  here.  Because  the  degree 
of  ionization  of  the  flow  is  very  small  and  frozen,  the  equations 
governing  the  flow  properties  arc  unaffected  by  the  presence 
of  the  ionized  particles  and  the  problem  can  be  divided  into 
two  parts,  the  hypersonic  flow  about  the  sharp  leading  edge 
plate  and  the  behavior  of  the  ionized  particles  in  this  flow. 

A.  Hypersonic  Flow  About  a  Sharp  Leading  Edge  Flat  Plate 

As  discussed  often  in  the  literature  (for  example.  Ref.  12  ), 
the  flow  field  over  a  sharp  leading  edge  body  traveling  »t 
hypersonic  speeds  may  be  separated  into  four  distinct  regions 
(Fig.  l).  The  first  region,  the  kinetic  flow  regime,  is 
located  in  the  immediate  vicinity  of  the  sharp  leading  edge. 

In  this  region,  where  property  transfer  is  conducted  at  a 
molecular  level,  one  must  use  the  Boltzmann  equation  to 
adequately  describe  the  flow  field.  This  region  gives  way  to 
a  small  transition  region  in  which  the  flow  field  is  becoming 
increasingly  denser  and  hence  more  dcscribablc  on  a  continuum 
basis.  Continuing  downstream,  one  enters  the  merged  layer. 

The  merged  layer  is  characterized  as  a  fully  viscous  region 
in  which  no  discernible  distinction  can  be  made  between  the 
boundary  layer  and  the  shock  layer.  Also  characteristic  of 
this  region  are  the  velocity  and  temperature  slips  that  occur 


nt  a  wall  and  a  t».nall  region  of  rarefied  flow  in  the  upper 
portions  of  the  viscous  layer. 

Progressing  further  downstream,  the  flow  asymptotes  into 
the  interaction  reqime.  The  initial  portion  of  this  region 
is  referred  to  as  a  strong  interaction  region  in  which  a 
snail  inviscid  layer  is  being  formed.  However,  there  exists 
a  strong  coupling  between  this  inviscid  layer  and  the  adjacent 
viscous  layer  such  that  one  cannot  treat  each  layer 
independently.  It  is  only  until  one  proceeds  further 
downstream  into  the  weak  interaction  regime  that  the  usual 
boundary  layer  techniques,  as  first  described  by  Prandtl,  are 
applicable. 

It  becomes  desirable,  therefore,  for  one  to  develop  a 
set  of  continuum  governing  equations  which  are  uniformly  valid 
from  the  wall  to  the  free stream  and  adequately  describes  the 
flow  field  from  the  merged  layer  downstream  into  the  weak 
interaction  region.  Such  an  analysis  was  first  carried  out 
by  Rudman  and  Rubin,  (Ref.  13)  in  which  they  postulated  one 
set  of  equations  which  were  valid  throughout  the  disturbance 
region. 

Rudman  and  Rubin  started  with  the  governing  Navier-Stokes 
equations  non-d  intensions  li  zed  with  respect  to  local  reference 
values.  They  postulated  that  if  the  shock  layer  on  a  slender 
body  is  considered  small,  such  that  the  normal  gradients  are 
much  larger  than  either  of  the  two  lateral  planar  gradients 
(this  ratio  being  denoted  as  *"*)  one  may  write  a  first  order 
approximation  to  the  full  Navier-Stokes  equations.  An 
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exception  occurred  in  the  energy  equations  where  6®  terms 

were  retained  in  order  to  obtain  a  non-trivial  solution  (Ref.  . 

Tref 

Simultaneously,  a  second  parameter  A?=  ' '^*"5  appeared  in  the 

00 

expanded  equations.  Since  Trgf  represents  the  order  of 
magnitude  of  the  local  temperature,  one  can  easily  show  that 
0.5  >  A2  >  .12  depending  upon  the  wall  temperature  (that  is, 
cold  or  adiabatic  wall)  and,  therefore,  may  also  be  considered 
small.  Consequently,  A®  terms  were  also  neglected  which  lead 
to  the  elimination  of  the  streamwise  pressure  gradient.  The 
consequences  of  this  are  discussed  in  Ref.  13  .  Within  this 
approximation  one  can  treat  the  x-momentum  equation  as  a 
parabolic  partial  differential  equation  in  which  the  inertia- 
viscous  terms  balance.  This  immediately  makes  possible  the 
use  of  a  finite  difference,  initial  value  technique  to  obtain 
the  flow  field  description. 

One  may  write  this  set  of  parabolic  equations  for  a  two- 


dimensional  flow  as  follows: 

Continuity:  {s*)}  (!) 

x-momentum :  />UUW*  Ug)|  (2) 

y-momentum:  +  -  -  Pj  (/Wg]  +(>  Ug]x  -  3*  (3) 

Energy.  —  +  (4) 

state:  p»^7  (5) 


where  p, p,T,  and  u  are  non-dimensional ized  with  respect  to 

their  respective  freest ream  values;  w  with  ;  v  with  u 

6  i  yM.su 

x  with  L  and  y  with  6.  Note  that  6=  =-  =  r-  -  -  and  L=  — — — 

L  y^M  p»u» 

CD 
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For  a  detailed  derivation  and  an  in-depth  discussion  of  these 

equations,  one  should  refer  to  Ref.  13  . 

The  following  initial  and  boundary  conditions  were  used 

in  the  numerical  solution  for  the  problem  in  question.  The 

initial  profiles  used  were  the  uniform  freestream  conditions. 

It  was  shown  by  Rudman  and  Rubin  that  their  solution  of  the 

flow  is  not  sensitive  to  the  initial  condition  for  V  <0.4. 

00 

For  larger  values  of  VM  the  solution  of  the  flow  depends  on 
the  initial  conditions  applied  at  the  leading  edge. 

Slip  boundary  conditions  are  enforced  on  the  surface  of 
the  plate;  namely, 

Y  =  O  X>  O 

U=AUj,  V=0  ,  T.  Tw  +  ^-  -£-T3  (6) 

Y  — ►  00  X  >  O 


V  — -  o 

A  numerical  solution  of  these  equations  was  carried  out 

for  the  flow  conditions  in  question  on  a  CDC  6600.  The 

conditions  are  as  follows:  *^=18.0,  Tq=7560oR,  T^T^O.075, 

*k 

a  =  0 ,737 ,  Y=1 .4,  c^=0.85. 
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B .  Charged  Particle  Density  Profiles 


In  general,  the  ionized  flow  about  a  body  can  be  divided 

into  an  inviscid  and  viscous  region;  the  viscous  layer  is  then 

divided  into  two,  the  sheath  and  ambipolar  region  (Ref. 14-16). 

In  the  latter,  convection  and  diffusion  operate  simultaneously, 

and  charge  neutrality  is  maintained.  The  thickness  of  this 

diffusion  layer  is  of  the  same  order  of  magnitude  as  the 

viscous  momentum  layer,  i.e.,  of  order  .  In  the  sheath, 

,/Re 

charge  separation  takes  place*  convection  is  unimportant  .  The 
main  transport  mechanisms  in  the  sheath  are  diffusion  and  mobil¬ 
ity  of  ions  and  electrons.  The  thickness  of  this  region  is 
about  10  times  the  Debye  shielding  distance,  evaluated  from 
plasma  properties  at  the  sheath  edge. 

Under  the  flow  and  plasma  conditions  about  the  plate  in 
this  work,  the  sheath  is  much  thinner  than  the  disturbed 
region,  which  includes  the  viscous  and  outer  compression 
regions;  thus  charged  particle  diffusion  across  the  layer  was 
assumed  to  be  essentially  ambipolar  in  nature. 

Following  the  formulation  by  Chung  (Ref.  16)  the  charge  specie 
conservation  equation  in  the  quasi-neutral  region  of  the 
disturbed  layer  about  the  plate  can  be  written  as 

In  deriving  this  equation,  charge  neutrality  N^=$Je  has  been 
invoked  and  diffusion  in  the  x-direction  neglected. 
Nondimensionalizing  this  equation  in  the  same  manner  as  the 
governing  equation  for  the  hypersonic  flow  about  the  plate 


C  C 


(8) 
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where  p,  T.,  T  ,  u  and  c  are  non-dimensionalized  with  respect 
1  e 


to  their  respective  freestream  value?  with 


u  Is! ?  x  with  L  and  y  with  f,  and  d= 

CO 

r 

>>uc„  .^iljc(u|-) 


p  D. 

^  CO 


M 


v  with 
,  we  obtain 


(9) 


Here,  S  =  is  the  ion  Schmidt  number.  The  assumed  boundary 

conditions  in  this  case  were 


>4=0  X  »  O  X  *  0 

- - (1Q) 

C  »  Cw  +  JkC3  C  —  1 

The  equation  governing  the  charge  particle  distribution  was 
solved  simultaneously  with  the  hypersonic  flow  about  the 
plate.  In  the  calculation  the  electron  temperature  Tg  was 
assumed  to  be  frozen  and  T^  the  ion  temperature  was  assumed 
equal  to  the  static  temperature  of  the  gas  at  any  point  in 
the  flow  field. 

The  choice  of  had  some  effect  on  the  charge  particle 
density  profile,  as  shown  in  Fig.  2  .  Cw=0  was  use<*  to 
agree  with  the  cylindrical  probe  measurements  of  the  ion  density 
near  the  plate  surface.  The  effect  of  the  ion  Schmidt 
number  (Sc^)  on  the  ion  density  distribution  in  the  disturbed 
region  about  the  plate  under  our  experimental  conditions  is 
exhibited  in  Fig.  3  .  The  ion  density  profile  appeared 
to  be  insensitive  to  the  value  of  Sc^. 
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III.  LANGMUIR  PROBES  IN  FLOWING  PLASMAS 


The  cylindrical  electrostatic  probe  has  been  found  to  be 
an  attractive  diagnostic  tool  for  low  density  plasma  flows  due 
to  its  excellent  spatial  resolution  and  relative  simplicity 
of  construction.  Recent  reviews  by  de  Leeuw  (Ref.  17)  and  Chen 
(Ref.  18)  gave  detailed  accounts  of  Langmuir  probe  theory  and  its 
application. 

One  of  the  most  realistic  collisionless  probe  theories 
published  to  date  is  that  of  Laframboise  (Ref.  4 )  who  has 
used  the  mathematical  framework  of  Bernstein  and  Rabinowitz 
(Ref.  3  )  to  calculate  entire  probe  characteristics  for  both 
cylindrical  and  spherical  probes  at  rest  in  a  two-temperature 
plasma.  Laframboise  assumes  that  the  plasma  consists  of  two 
species  of  charged  particles,  one  positive  and  one  negative, 
each  having  a  Maxwellian  velocity  distribution  with  its  own 
characteristic  temperature  and  that  whenever  a  particle  of 
either  species  strikes  the  probe  surface,  it  is  annihilated 
and  removed  from  the  problem.  Laframboise  calculations  are 
based  on  the  Vlasov  equations  for  two  species,  coupled  by 
Poisson's  equation  for  the  electric  field,  and  the  results  are, 
therefore,  not  expected  to  be  valid  unless  binary  collisions 
can  be  ignored  and  the  electric  field  consistent  with  Poisson's 
equation  is  the  dominant  influence  on  the  motion  of  the  charged 
particles. 

For  a  plasma  where  the  positive  particles  are  singly 
charged  ions  and  the  negative  ones  are  electrons,  Laframboise' s 
results  may  be  summarized  as  follows: 


For  x  ' 
P 

Ji 


-J 


Je  *  eN  ^ 

dm* 

x  >  0 

P 

Ji  -eN 

Je  «eN 

pH 

|2JTMe 

(U) 

(12) 


e-\  *f 

**  (*P,  RpAd,^“) 


(13) 

(14) 


When  the  probe  is  at  plasma  potential,  x  =0,  either 

P 

species  arrives  at  the  probe  surface  at  a  rate  given  by  its 

random  kinetic  flux,  ,  and  when  a  species  is  repelled,  its 

flux  to  the  surface  is  reduced  from  this  value  by  the 

Boltzmann  factor.  When  it  is  attracted,  on  the  other  hand, 

the  situation  becomes  more  complex.  When  the  electrons  are 

attracted,  their  flux  increases  over  the  random  kinetic 

value  by  the  factor  a  ,  calculated  by  Laframboise.  When 

*e 

the  probe  is  made  negative  so  as  to  attract  ions,  the  ion 

yux  soon  reaches  a  magnitude  of  the  Cider  eNVj- 5P  •  wh*ch  i* 

jjr’  times  the  random  ion  flux  and  may  be  much  larger  than 

the  latter  if  T  »T..  Again,  a  represents  a  dimensionless 
e  i  p^ 

correction  factor  which  emerges  from  theory. 

Both  r*  and  a  depend  not  only  on  the  dimensionless 
pi  pe 

potential  difference  between  the  probe  and  the  plasma,  x  , 

P 

but  also  on  the  properties  of  the  plasma  through  R  /X_  and 

P  13 

TiAe#  as  well  as  on  probe  geometry.  Both  are  generally  of 

order  unity  at  moderate  probe  potentials  and  values. 

The  magnitude  of  a  and  ?  for  cylindrical  and  spherical 

pi  ^e 


12 


probes  are  given  in  Ref.  4  for 

|Xpl  <  25 

0<  yTe  <  1  (15) 

o  <  <  101 

Figure  4,  which  is  reproduced  here  from  Laframboise  work, shows 

the  dependence  of  ~  on  R  / ^  for  T^/Te=0.  Included  in 

Pi  ?  i  e 

Fig.  4  is  the  orbital  motion  limited  current  boundary. 

Current  collection  is  said  to  be  orbital  motion  limited  at 


the  energy  E,  when  no  absorption  radius  exists  for  the  energy 
E  and  the  current  is  limited  by  the  large  angular  momentum  of 
the  particles.  The  orbital  motion  limited  current  represents 
the  maximum  current  of  particles  of  a  given  energy  that  can 
be  collected  in  the  collisionless  case  by  a  given  potential. 
From  this  figure,  it  is  seen  that  as  R  /\~  decreases,  the 

P  u 

current  c,  reaches  the  orbital  motion  limit  and  levels  off 
to  a  constant  value. 

Since  the  publication  of  Laframboise  results,  many 
experimental  investigations  were  carried  out  to  verify  the 
applicability  of  the  static  probe  theory  to  probes  in 
collisionless  flowing  plasmas  (Ref.  19)  .  In  his  work,  Sonin 
examined  the  behavior  of  free  molecule  cylindrical  Langmuir 
probes  in  supersonic  flows.  Those  experiments  indicated  that 
under  flow  conditions,  the  theory  of  Laframboise  represents 
the  ion  current  characteristics  adequately  at  large  R  /\D» 
Departure  from  theoretical  ion-currents  were  observed  in  the 
small  probe  limit,  when  the  ratio  of  cylindrical  probe  radius 
Rp  to  plasma  Debye  length  ^  was  less  than  about  3.  Similar 
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departure  was  observed  by  Lederman  et  al.in  their  experiments 

with  cylindrical  probes  in  rarefied  ionized  hypersonic  flow 

(Ref.  20)  .  The  leveling  off  of  the  a  curves  for  small  values 

^i 

R  /X  in  La  framboise  results  (see  Fig.  4  )  was  not  compatible 

p  D 

with  those  experimental  results,  a  ,  exceeded  markedly 


S' 


the  theoretical  collisionless  values  at  low  R  /XD. 

In  many  experiments  by  Lederman  e_  al.  conducted  in  the 
hypersonic  shock  tunnel,  the  parameter  R^/X^  was  much  smaller 
than  unity.  This  resulted  from  the  combined  effect  of  the 
plasma  conditions  (N.  ,Te)  and  small  probe  radii  used  to 

CO 

minimize  gasdynamical  effects  on  collisions.  Also  in  the 
course  of  this  work,  length  limitations  required  for  the 
achievement  of  spatial  resolution  in  ion  density  distribution 
measurements  resulted  in  length -to-diameter  ratio  //2Rp  <  100. 
It  appeared  that  experimental  information  on  the  effect  on  ion 
current  collection  of  the  parameter  £/2Rp<100  and  Rp/XD<  3  in 
rarefied  hypervelocity  ionized  flow  is  needed. 

This  lead  to  an  experimental  investigation  by  Lederman 
et  al. (Ref.  5)  which  determined  the  behavior  of  biased 
electrostatic  probes  in  the  regime  of  high  rarefaction,  and  low 
degree  of  ionization  hypervelocity  flow.  Using  the  slightly 
ionized  flow  in  a  hypersonic  shock  tunnel,  data  for  various 
diameter  and  length-to-diameter  ratio  ( «t/2R^ probes  was  obtained 

under  known  plasma  conditions  (N.~  3  x  10eel/fcm3,  T  »  3000°K, 

1®  e_ 

CO 

T^/Te~U) .  The  free  stream  ion  number  density  was  measured 
with  microwave  resonant  cavities  and  the  electron  temperature  was 
obtained  using  cylindrical  probes  operating  in  the  electron 
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retarding  field  region.  These  experimental  results  showed 
that  the  current  collection  ability  of  a  probe  operating  in 
the  above-mentioned  regime  and  the  resulting  current  density 
is  greatly  enhanced  with  the  decrease  of  the  diameter  of 
the  cylindrical  collecting  probe  and  the  subsequent  decrease 
of  the  ratio.  This  effect  thus  presents  a  significant 

departure  from  Laframboise  theory  at  the  point  where, 
according  to  that  theory,  the  current  should  reach  its  orbital 
motion  limited  value  and  should  no  longer  increase.  The 
experimental  values  of  the  normalized  nondimensional  current 
a  for  R  /XD  as  obtained  by  Lederman  et  al.  are  presented  in 
Fig.  5  .  The  observed  experimental  increase  in  the  ion  current 
at  small  value  of  R^/\D  can  be  attributed  to  collisions  of 
the  orbiting  ions  as  pointed  out  by  Shih  and  Levi  (Ref.  21). 
Since  the  Coulomb  cross  section  for  ion-ion  collisions 
increases  rapidly  as  the  ion  temperature  decreases,  in  the 
limit  of  small  ion  temperature,  which  was  the  case  in 
Lederman' s  work,  there  always  exists  ion-ion  collisions  in 
sufficient  numbers  to  destroy  orbital  motion.  Shih  showed 
in  this  work  that  there  is  an  increase  in  current  collection 
by  the  probes  due  to  small  collisional  effects  in  the  region 
where  orbit  motion  limit  should  have  been  reached  in  a 
collisionless  plasma  as  predicted  by  Laframboise. 

Lederman  et  aL  also  formed  a  relation  between  current 
density  and  the  length-to-diameter  ratio  as  shown  in  Fig.  6  , 
thus  permitting  the  use  of  very  small  probes  to  achieve  better 
spatial  resolution.  The  increase  in  current  density  collected 
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by  probes  with  //2R  <100  can  be  attributed  to  end  effects  when 
the  Debye  length  is  large  compared  with  the  probe  radius  and 
the  flow  speed  is  high.  Under  those  conditions,  the  ion 
current  collected  by  the  probe  consists  of  ions  entering  the 
side  sheath  arid  of  ions  entering  through  the  front  of  the 
sheath.  For  short  enough  probes,  the  latter  dominates,  thus 
causing  the  large  increase  in  observed  current  collected  by 
short  probes. 
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IV.  ELECTRON  BEAM  DENSITY  PROBE 


The  optical  electron  beam  density  probe  was  suggested 
and  developed  by  Schumacher  and  Gadamer  (Ref.  22)  and  was 
first  applied  to  an  actual  flow  investigation  by  Gadamer 
(Ref.  23) #  who  showed  that  it  can  be  used  to  accurately 
determine  the  local  densities  in  a  rarefied  gas  flow.  The 
technique  employs  a  narrow,  well-collimated,  high  energy, 
electron  beam  which  is  projected  at  right  angles  to  the  flow 
through  the  region  of  interest.  Because  of  the  relatively 
small  collision  cross-section  of  the  gas  atoms  for  20  keV 
electrons,  the  interaction  of  the  gas  with  the  beam  is  weak 
and  results  mainly  in  the  scattering  at  small  angles  of  some 
of  the  primary  electrons.  At  the  electron  energy  used,  many 
of  the  electron-atom  interactions  are  of  the  elastic  type. 

There  occur,  however,  sufficient  inelastic  collisions, 
resulting  in  electronic  excitation  and  subsequently  de¬ 
excitation,  to  give  rise  to  luminescence,  which  is  largely 
confined  to  the  beam  region  itself.  By  focusing  the 
luminescence  onto  an  optical  slit,  which  is  at  right  angles 
to  the  beam,  the  light  issuing  from  a  small  elemental  gas 
volume  is  selected  and  can  subsequently  be  measured  through 
the  utilization  of  a  photomultiplier.  Under  certain  conditions, 
the  photor>current  measured  is  a  direct  indication  of  the 
molecular  number  density  in  the  elemental  volume  under 
consideration. 

It  should  be  pointed  out  that,  at  gas  densities  corresponding 
to  a  gas  pressure  of  200uHg  at  room  temperature,  the  average 
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time  between  atomic  collisions  of  order  10  sec.  for  most 
gases  becomes  comparable  to  the  lifetimes  of  excited  states 
of  the  atoms  of  order  10  sec.  for  allowed  transitions. 

Under  these  conditions,  a  measurable  fraction  of  the  electron 
beam  excited  atoms  collide  with  other  atoms  before  they  can 
emit  the  excitation  energy  in  the  form  of  radiation.  This 
makes  the  observed  intensity  of  the  electron  beam  luminescence 
less  than  one  would  otherwise  predict.  The  effect  is  generally 
referred  to  as  collisional  quenching.  Since  the  collision 
frequency  and,  therefore,  the  collisional  quenching  effects 
are  temperature  dependent  as  well  as  density  dependent,  it  is 
almost  imperative  to  work  with  gas  densities,  at  which 
collisional  quenching  can  be  neglected.  For  the  visible 
spectrum  of  nitrogen,  Gadamer  (Ref. 23)  found  no  quenching 
effects  up  to  lOOuPg. 

In  addition,  there  exists  the  possibility  of  the 
luminescence  being  carried  downstream  in  a  high  velocity  gas 
flow.  As  pointed  out  before,  for  the  nitrogen  first  negative 

.7 

system  the  lifetime  is  approximately  10  sec.  Since  the 
velocity  in  the  HST  is  approximately  3  x  10s  cm/sec.,  the 
drift  in  this  case  is  about  0.03  cm,  and  it  is,  therefore, 
possible  to  obtain  good  spatial  resolution  by  observing  the 
center  portion  of  the  beam  having  a  2 -mm  diameter. 

A.  Origin  of  Luminescence  in  Nitrogen 

The  electron  beam  technique  involves  the  ionization  and 
excitation  of  nitrogen  molecules  with  a  beam  of  highly 
energetic  electrons  (10-30  keV) .  The  resulting  spontaneous 
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emission  in  the  spectral  region  3500-5000  A  is  predominantly 
that  of  the  first  negative  system  of  nitrogen  (1-) . 

Analysis  by  Muntz  (Ref.  24  and  Petrie  (Ref.  25)  proposed 
the  excitation  process  to  be  directly  from  the  neutral  ground 
state  molecules  (NjX1?*)  to  the  excited  states  of  the  molecular 
ion  (N^B*^)  by  inelastic  collisions  with  high  energy  electrons 
(See  Fig.  7).  This  analysis  discounted  other  paths  of 
excitation  to  such  as  cascade  population,  double 

excitation,  and  secondary  excitation  of  ground  state  ions. 

The  observed  emission  corresponds  to  the  electronic 
transition  M^B*^  -.  (the  ground  state  of  the  molecular 

ion)  with  the  average  lifetime  of  the  excited  ions  on  the 
order  of  10“7  seconds.  Muntz  (Ref.  24)  also  showed  that  for 
low  density  gas  flows  at  room  temperature  the  excited  gas 
particles  should  experience  little  interference  due  to  gas 
kinetic  collisions  during  the  process  of  excitation  and 
emission. 

The  luminescence  caused  by  the  de-excitation  of  the 
molecules  can  be  examined  in  the  following  manner.  The  rate 
of  emission  of  photons  of  frequency  ^  corresponding  to 
radiative  transition  from  molecular  state  E ^  to 
of  the  luminescing  gas  is 

$(V* k)  •  A)k  Nj  (16) 


per  cm 


where  Nj  is  the  density  of  molecules  per  cm3  in  the  parent 
state  Ej  and  A^  is  the  spontaneous -emission  coefficient  for 
the  transition.  Under  the  weak  excitation  condition  employed 
here,  is  produced  by  binary  encounters  between  beam 
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electrons  and  neutral  ground  state  Nj  moleculesi  the  rate  of 
production  of  N.  is  then 

3  (j?)wc‘«»(£)pl('r)  (17) 

where  E  is  the  energy  of  the  incident  electrons, P  the  pressure 

of  Nj#  i(r)  the  spatially  dependent  flux  of  beam  electrons, 

and  k.(E)  is  related  to  the  cross-section  for  excitation  of 
3 

the  Ba  v  state  of  by 

Kj  (e)  =  2.22  *  1029  (e)  (18) 

The  Ba£  state  is  de-excited  in  two  ways.  First,  the  molecule 

can  radiate  to  all  lower  states  yielding  a  spectrum  which 

o 

includes  the  3914  A  line  (the  most  intense  radiation)  and  the 


rate  of  radiation  is  given  by 


m 


0-1 

*  -I  Ajh  Nj  =  *  AjNi 
RAP  K*0 


(19) 


where  A  ^  is  the  reciprocal  of  the  radiation  lifetime  of  state 

N.  and  the  radiation  lifetime  of  state  N.  is  the  average  time 
3  3 

spent  by  the  molecules  in  state  Nj . 

The  excited  molecule  can  also  lose  its  energy  in 
collisions  with  other  molecules*  this  process,  known  as 
quenching, is  described  by 

'CnPNj  (20) 


\dt  /  QUtN. 


Here  Cnis  the  quenching  coefficient  for  nitrogen.  In  the 
steady  state,  the  production  of  Nj  is  equal  to  its  total 


removal  rate 


hence. 


NjA)  «■  C„PNj  *  Kj  (e)PL  (?) 
*(p)  •  A>*Nl  ■  (r) 


(21) 


(22) 
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Since  for  a  constant  temperature,  P  is  prcpor ' ional  to  the 
number  of  Nj  particles  present,  Eq.  (22)  shows  that  there 
exists  a  linear  relation  between  the  radiant  intensity  and 
the  density  of  nitrogen  as  long  as  quenching  effects  are 
small  (Cn«  Aj).  This  linear  relationship  is  the  basic 
principle  of  the  electron  beam  density  probe. 
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V.  EXPERIMENTAL  APPARATUS 


A.  Experimental  Facility 

The  experiments  reported  throughout  this  work  were 
performed  in  the  hypersonic  shock  tunnel  of  the  Polytechnic 
Institute  of  Brooklyn#  located  at  the  Graduate  Center  in 
Farmingdale.  The  schematic  layout  of  the  H.S.T.#  together 
with  its  principal  components,  is  shown  in  Fig.  8  .  The  HST 
employs  a  variable-area  shock  tube  where  the  driver  is  6  inches 
I.D.  and  24  feet  long,  and  the  driven  is  3  inches  I.D.  and 
28  feet  long.  The  driver  gas  was  helium  at  room  temperature 
and  the  driven  gas,  air  or  nitrogen.  The  shock  processed  gas 
at  the  end  of  the  shock  tube  is  expanded  through  a  conical 
nozzle  into  the  test  section.  A  uniquely  designed  valve 
located  at  the  end  wall  of  the  shock  tube  prevents  the  driver 
gas  from  flowing  into  the  test  section,  thus  avoiding  damage 
and  contamination  of  the  instrumentation.  A  detailed  description 
of  the  facility  and  its  performance  capability  is  given  in 
Ref.  26  . 

B.  Flow  Conditions 

The  operating  initial  conditions  of  the  shock  tunnel  were 
fixed  throughout  this  investigation.  With  an  initial  driven 
pressure  of  38  mmHg  air  or  nitrogren  and  driver  pressure  of 
1800  psia  helium,  a  shock  Mach  number  of  about  7  was  obtained. 
This  resulted  in  a  pressure  of  about  350  psia  and  temperature 
of  4000°K  behind  the  reflected  shock.  The  shock  processed 
gas  was  then  expanded  through  a  10°  half-angle  and  52-inch 
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exit  diameter  conical  nozzle  into  the  test  section.  The  flow 
conditions  in  the  test  section  were  obtained  in  the  following 
way.  The  stagnation  conditions  behind  the  reflected  shock 
were  obtained  from  shock  tube  tables  for  equilibrium  air  or 
nitrogen  for  the  measured  shock  velocity,  and  the  flow 
properties  in  the  nozzle  and  the  test  section  were  calculated 
using  the  CAL  program  for  nonequilibrium  nozzle  flows  (Ref.  27). 
Calculation  of  the  flow  properties  using  this  program  showed 
that  the  flow  becomes  frozen  at  about  25  in. downstream  from 
the  . nroat  for  stagnation  conditions  To«4000°K  and  P^SSO  psi. 

Table  I  lists  the  flow  conditions  at  positions  along  the 
nozzle  and  Table  II  lists  those  in  the  test  section.  At  the 
various  stations  along  the  nozzle  a  rake  of  cylindrical 
electrostatic  probes  aligned  parallel  to  the  flow  direction 
was  used  to  determine  the  uniformity  of  the  flow.  Stagnation 
heat  transfer  and  ion  density  surveys  made  under  test 
conditions  in  the  test  section  showed  the  flow  is  uniform 
through  a  core  26  inches  in  diameter.  This  region  is  only  a 
portion  of  the  nozzle  exit  diameter  because  of  boundary  layer 
growth  in  the  nozzle  due  to  the  low  Reynolds  number  flow.  The 
relatively  large  uniform  region  allows  one  to  introduce  fairly 
large  models  for  surface  measurements  and  flow  field  investigation. 

C.  Electron  Beam  Density  Probe 

The  electron  beam  density  probe  technique  is  based 
on  the  stimulation  of  emission  from  a  gas  by  collisional 
excitation  using  a  beam  of  energetic  electrons  (larger  than 
10  keV) .  Inelastic  collisions  between  the  gas  molecules  and 
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the  beam  electrons  result  in  electronic  excitation  and  sub¬ 
sequently  de-excitation  giving  rise  to  luminescence  which  is 
mostly  confined  to  the  beam  region.  The  local  emitted 
intensity  is  directly  proportional  to  the  local  gas  number 
density  of  molecules.  With  the  luminescent  light  yield  detected 
at  different  points  along  the  beam  with  a  photomultiplier, 
spatial  density  measurements  can  be  carried  out.  An  electron 
density  probe  basically  consists  of  two  systems;  one.  the 
electron  gun  which  is  used  to  generate  a  narrow  collimated, 
high  energy  electron  beam. and  the  other  is  the  electro-optical 
system  for  point-wise  luminescence  detection.  The  electron 
density  probe  used  in  this  investigation  is  described  in  the 
following  section. 

D.  Electron  Gun  and  Electro-Optical  Systems 

The  electron  gun  assembly,  complete  with  diffusion  pump 
and  gauges,  is  mounted  on  the  hypersonic  shock  tunnel  as  shown 
in  Fig.  9  .  An  external  view  of  the  gun,  battery,  power 
supplies,  and  instrumentation  is  shown  in  Figs.  10  and  11. 

The  electron  beam  is  produced  by  a  specially  designed 
electron  gun  ,  which  consists  of  a  tungsten  filament  and  one 
accelerating  grid  (the  filament  is  a  plug-in  type  allowing 
easy  replacement) .  The  advantage  of  a  tungsten  filament 
electron  gun  over  a  commercial  electron  gun  with  a  Barium 
Cathode  is  twofold.  First,  the  filament  can  be  exposed 
to  atmospheric  pressure  without  losing  its  emission  capability, 
and  second,  no  differential  pumping  is  needed  to  maintain  high 
vacuum  in  the  gun  chamber  since  it  can  operate  at  higher 
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pressure,  about  10~4  Torr,  compared  to  at  least  10~  required 
for  a  commercial  gun.  The  gun  is  housed  in  a  chamber  constantly 
evacuated  by  a  230  liter/sec.  oil  vapor  diffusion  pump  (Kinney 
KDB-200)  and  maintained  at  approximately  10  Torr  pressure. 

The  pressure  in  the  gun  chamber  is  monitored  by  an  ionization 
gauge  (Veeco  RG-31A) .  In  order  to  reduce  oil  deposits  in  the 
gun  chamber  caused  by  backstreaming,  the  diffusion  pump  was 
equipped  with  a  nitrogen  cold  trap  (Ki.iney  KDB-200).  A  special 
diffusion  pump  oil  with  low  backstreaming  properties  was 
used  (Convslex-10  manufactured  by  Bendix) . 

From  the  gun  chamber,  the  electron  beam  passes  into  a 
1-inch  diameter  ard  an  18-inch  long  drift  tube,  and  from  there 
into  the  test  section  through  an  aperture  of  low  gas 
conductance.  This  restricts  the  inflow  of  gas  from  the  test 
section  and  confines  the  electron  beam  geometrically.  The 
aperture  consists  of  a  cylindrical  hold  2  mm  in  diameter  and 
1  inch  long  drilled  into  a  brass  insert,  which  can  easily  be 
replaced. 

Careful  alignment  and  focusing  of  the  beam  is  necessary 
to  allow  it  to  pass  through  the  aperture  into  the  test  section. 
This  is  achieved  by  an  electromagnetic  focusing  and  deflection 
system.  Fig.  12.  This  system,  which  consists  of  two 
deflection  coils  and  one  focusing  coil  (Celco  AF  304-200)  is 
mounted  on  the  drift  tube.  Regulated  power  supplies  (KEPCO) 
were  used  to  supply  the  control  current  to  the  coils.  After 
passing  through  the  test  section  the  electron  beam  impinges 
on  a  brass  block  which  serves  as  the  beam  collector.  The  beam 
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collector  is  electrically  isolated  from  ground  by  a  teflon 
insulator,  so  that  the  collected  beam  current  can  be  measured. 
The  current  to  the  filament  is  supplied  by  a  6  volt  heavy 
duty  battery  and  can  be  adjusted  by  a  rheostat  to  control 
the  emission.  The  battery  which  is  at  beam  potential  is 
placed  in  a  big  plexiglass  box  because  of  the  high  voltage 
with  respect  to  ground. 

The  cathode  potential  is  supplied  by  a  Universal  Voltronic 
high  voltage,  unregulated,  D.C.  power  supply  which  can  be 
varied  between  zero  and  32  KV.  The  anode  is  kept  at  ground 
potential  or  a  few  volts  positive  bias  with  respect  to  ground. 
The  grid  voltage  is  obtained  from  a  voltage  divider  circuit, 
which  allows  ample  adjustment  of  the  beam  current  (see 
Fig.  12)  . 

The  function  of  the  electro-optical  system  is  to  measure 
the  intensity  of  the  radiation  issuing  from  a  small  gas  volume 
in  the  electron  beam.  In  accordance  with  this  requirement, 
the  instrument  must  have  a  high  sensitivity  to  detect  a  very 
weak  light  source,  and  it  must  spectrally  discriminate 
between  the  light  emitted  by  the  nitrogen  due  to  excitation 
by  the  beam  electrons  and  other  sources  of  luminescence 
present  during  a  test  in  the  hypersonic  shock  tunnel. 
Furthermore,  only  light  which  originates  in  a  well-defined 
region  of  the  electron  beam  should  be  allowed  to  be  picked  up 
by  the  photomultiplier. 

Good  sensitivity  is  obtained  by  using  an  optical  system 
of  low  f-number  in  conjunction  with  a  high  gain  photomultiplier. 
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Two  achromatic  lenses,  £  3.7  and  f  5.8,  are  used  to  form  a 

real  image  of  the  electron  beam.  In  the  image  plane  an  iris 

is  used  to  spatially  define  the  section  at  the  electron  beam, 

from  which  light  is  received.  A  narrow  band-pass  filter  is 

o 

used  to  select  the  most  intense  nitrogen  line  (3914A)  and  the 

intensity  of  the  transmitted  light  is  then  measured  with  an 

uncooled  photomultiplier  EMI  9502S  having  an  S-ll  spectral 

o 

response  (a  peak  at  400QA) .  The  spectral  response  curve  of 
the  photomultiplier  and  the  transmission  curve  of  the  narrow 
band-pass  filter  are  given  in  Fig.  13. 

It  should  be  noted  that  the  lens-to-beam  and  lens-to-iris 
distances  are  approximately  192  and  275  mm,  respectively,  giving 
magnification  close  to  1.5.  Thus  the  3  mm  diameter  iris 
opening,  when  imaged  on  the  2  mm  wide  electron  beam,  is  only 
2  mm  in  diameter.  This  defines  the  elemental  volume  in  the 
beam  luminescence  where  the  measurements  were  made. 

The  photomultiplier  tube  and  the  optics  are  housed  in 
a  container  made  of  aluminum  tubing  (Fig.  14).  This  combined 
electro-optical  system  is  mounted  on  a  frame  inside  the  test 
section  (see  Fig. 9  ),  which  by  itself  is  mounted  rigidly  to 
the  electron  gun  housing.  This  arrangement  has  the  advantage 
that  optics  remain  aligned  with  the  electron  beam.  The 
photomultiplier  tube  and  the  optics  are  moved  from  point-to- 
point  along  the  beam  by  a  remotely  controlled  driving  system 
and  the  position  of  the  system  is  monitored  by  a  potentiometer- 
ammeter  readout.  One  of  the  problems  encountered  by  che  fact 
that  the  photomultiplier  is  inside  the  test  section  where  high 
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vacuum  is  maintained,  was  electrical  breakdown  between  the 
high  voltage  connections  leading  to  the  photomultiplier  at 
about  50  «Hg.  Suitable  precautions  have  been  taken  to  prevent 
this  from  happening. 

The  cathode  voltage  to  the  EMI  9502S  photomultiplier  was 
supplied  by  a  highly  regulated  D.C.  power  supply  (N.J.E.  S-326) . 
Sec  Fig.  14  for  the  tube  circuitry.  The  anode  current  was 
displayed  on  a  Tektronix  565  oscilloscope  and  recorded  by  a 
Polaroid  camera.  This  kind  of  recording  system  is  used  because 
of  the  very  short  test  duration. 

In  order  to  obtain  simultaneously  density  measurements 
at  two  points  during  one  test,  a  second  electro-optical  system 
almost  identical  to  the  first  one,  was  constructed.  Both 
systems  were  arranged  in  such  a  way  that  their  focal  points 
intersect  with  the  beam  and  the  distance  between  the  focal 
points  was  kept  at  1"  during  all  tests  of  the  near  wake  density 
survey,  and  0.2"  for  all  the  tests  in  the  flew  field  about  the 
plate.  Both  systems  were  connected  rigidly  to  a  plate  which 
could  be  moved  along  the  beam  with  the  driving  system  designed 
for  this  purpose. 

E.  Langmuir  Probes  and  Associated  Circuitry 

Cylindrical  Langmuir  probes  were  made  by  allowing  a 
length  of  tungsten  wire  to  protrude  from  a  piece  of  Teflon 
tubing  used  as  insulation,  having  an  inside  diameter  of 
approximately  the  diameter  of  the  wire.  Tungsten  wires  0.001- 
0.080  inches  in  diameter  were  used  for  all  probes.  With  the 
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desired  length  of  wire  protruding  from  the  tube,  the  interface 
between  the  Teflon  tube  and  the  tungsten  wire  was  filled  with 
an  epoxy,  thus  making  a  firm  bond.  This  process  eliminated 
any  possibility  of  faulty  current  leakage  through  the  Teflon- 
wire  interface.  A  stainless  steel  tube  was  used  as  a  protective 
outer  shield  for  the  probe.  To  facilitate  Langmuir  probe 
boundary  layer  surveys  over  the  entire  flat  plate  surface, 
a  group  of  25  identical  probes  0.020  inches  in  diameter 
and  0.5  inches  long  was  arranged  in  a  rake  for  simultaneous 
acquisition  of  data.  The  probes  were  of  the  "plug-in"  type 
allowing  easy  replacement  or  interchanging. 

A  second  rake  of  probes  was  constructed  for  the 
investigation  of  the  effect  of  density  increase  on  ion  current 
collected  by  cylindrical  electrostatic  probes.  The  rake 
consisted  of  four  rows  of  eight  probes  each.  One  row  contained 
eight  probes  of  the  same  diameter  and  length  to  check  the 
uniformity  of  the  flow  and  the  next  two  rows  contained  two 
sets  of  probes  of  the  same  -t/2R  ratio  but  of  variable  diameter 

Xr 

ranging  from  0.001  inch  to  0.080  inch.  The  last  row  consisted 
of  probes  with  same  diameter  but  different  length  to 
investigate  fineness  ratio  effects. 

Lang».uj.r  probes  were  cleaned  ultra sonically  ii  acetone 
solution  before  being  installed  in  the  rake.  Probe  current 
collection  appeared  to  be  very  reproducible  in  different 
tests  for  the  same  test  conditions,  thus  discounting  any 
possibility  of  contamination  after  installation  and  between 
tests. 
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Figure  1  f>  allows  the  construction  and  circuit  for  the 
I.ur.qmmr  probe.  Prole  voltage  was  supplied  by  a  regulated 
D.C.  power  supply.  During  all  the  tests  reported  here, 
probes  were  biased  at  -9  volts  tor  ion  collection.  The  probe 
current  dm  inq  a  test  is  displayed  on  a  565  Tektronic  oscilloscope 
and  recorded  with  a  Polaroid  camera. 

I' .  Flat  Plates  and  Imbecded  Flush  Electrostatic  Probes 

Two  completely  different  plates  were  used  in  this  work. 

The  first  was  made  of  aluminum  and  a  stainless  steel  sharp 
leading  edge  (0.005")  as  shown  in  Fig.  16.  This  plate  was 
used  for  the  investigation  of  the  slightly  ionized  flow  in 
the  disturbed  region  over  the  plate.  Holes  were  bored  into 
the  surface  to  accommodate  the  flush  probes.  Four  rows  of 
those  probes  were  imbedded  flush  with  the  surface.  The  probes 
were  circular  in  shape,  all  made  out  of  brass,  insulated  from 
the  plate  with  plexiglass  sleeves,  and  varied  in  size  from 
1/16"  to  3/4"  diameter.  After  installing  the  probes  in  the 
plate,  the  combination  plate  probes  were  then  ground  for  a 
smooth-surface  finish.  Again,  the  flush  probes  were  cleaned 
before  installing  the  plate  in  the  test  section,  and 
reproduction  between  tests  was  very  good.  The  plate  was 
mounted  on  a  movable  string  to  allow  one  to  position  the 
plate  at  different  axial  locations  with  respect  to  the  stationary 
cylindrical  probe  rake. 

The  second  plate  was  used  for  flow  field  investigation 
with  the  electron  beam  density  probe.  Figs.  17  and  18.  Actually* 
two  plates  were  used  for  this  purpose;  the  first, made  of  1/16  inch 
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aluminum  sheet  with  a  15°  bevel  on  the  undersides, was  12  inches 
by  18  inches.  This  plate  was  used  to  obtain  density  profiles 
at  the  same  axial  positions  where  ion  density  profiles  were 
already  measured  with  the  electrostatic  probe.  The  second 
plate  was  made  of  0.015  eh  thick  spring  steel  sheet  6 
inches  by  6  inches.  This  plate  was  used  for  density 
measurements  in  the  vicinity  of  the  leading  edge  (x=2.0  inches). 
This  plate  was  introduced  for  measurements  near  the  leading 
edge  in  order  to  avoid  reported  observed  influence  of  a  finite 
thickness  plate  on  the  flow  field  development  near  the  leading 
edge  of  the  plate  (Ref.  28  ) .  These  plates  were  connected 
to  the  drift  tube  of  the  electron  gun  system,  and  a  2  mm 
diameter  orifice  machined  into  the  plate  allowed  the  beam  to 
emerge  through  the  plate  and  impinge  on  the  beam  collector. 

Since  the  electron  beam  density  probe  was  fixed  in  position 
in  the  tunnel,  measurements  were  made  starting  with  the  rear 
stations  and  moving  forward.  As  each  station  was  completed, 
a  new  orifice  was  machined  at  the  next  location  forward.  This 
procedure  prevented  any  disturbance  to  the  measuring  region 
from  orifices  ahead  of  the  testing  location. 

G.  The  Cone 

As  mentioned  previously,  density  distribution  in  the 
near  wake  of  a  10°  half  angle  4"  base  di.  meter  cone  was  also 
investigated.  The  cone  was  made  of  aluminum  with  a  stainless 
steel  tip  and  suspended  from  a  frame  with  0.020"  diameter  piano 
wires,  as  shown  in  Pigs. 19,20.  it  was  found  in  previous  work  of 
Lederman  (Ref.  29)  that  this  size  wire  docs  not  introduce  any 
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disturbance  in  the  wake  region.  The  frame  on  which  the  cone 
was  suspended  was  mounted  on  a  movable  string  to  allow  one  to 
position  the  cone  at  different  axiaj.  positions  with  respect 
to  the  electron  beam  density  probe. 
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VI.  EXPERIMENTAL  RESULTS 


A.  Test  Duration 

One  of  the  problems  generally  encountered  working  in 
hypersonic  shock  tunnels  is  the  test  juration  and  where  is 
the  test  in  the  time  history  of  the  recorded  data.  In 
previous  work  by  Lederman  et  al.  (Ref.  20)  and  in  the  present 
work,  the  test  duration  of  the  ionized  flow  in  the  hypersonic 
shock  tunnel  of  PIB  was  established  usina  cylindrical 
electrostatic  probes  and  microwave  cavities. 

Using  these  methods,  the  ionized  flow  under  the  present 
test  conditions  was  found  to  be  about  0.4  msec,  in  the  course 
of  the  present  investigation  using  the  electron  beam  density 
probe  and  stagnation  pressure  measurements,  a  uniform  flow 
region  of  a  longer  duration  which  starts  after  the  ionized 
flow  in  the  time  history  of  the  recorded  data  was  utilized. 
Those  two  flow  regimes  are  apparent  from  the  following 
figures  (Fig.  21  through  Fig.  23). 

In  Fig.2la,  two  traces  of  ion  current  collected  by 
cylindrical  electrostatic  probes  are  shown.  It  is  evident 
from  these  traces  that  after  the  initial  flew  containing 
ionized  particles  a  region  of  non-ionized  or  at  most  very 
slightly  ionized  flow  exists. 

In  Fig.  21b,  the  corresponding  background  radiation  of 
the  freestream  flow  without  beam  excitation  is  shown.  It  is 
interesting  to  note  that  radiation  occurs  at  the  exact  time 
of  the  ionized  flow  and  vanishes  simultaneously  with  the 
vanishing  of  the  ionization  of  the  flow. 
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In  the  next  figure.  Fig.  22,  two  traces  are  indicated. 

As  seen,  the  traces  of  both  oscillographic  recordings 
correspond  to  the  photoluminescence  in  the  freestream. 

(a)  corresponds  to  the  freestream  "off  beam"  luminescence 
whereas  (b)  corresponds  to  the  freestream  "on  beam"  luminescence . 
It  is  seen  that  the  initial  luminous  flux  corresponding  to 
the  ionized  portion  of  the  flow  is  followed  in  (a)  by  a 
non-radiating  flow  region  as  in  Fig.  21a  obtained  under 
"off  beam"  conditions  and  in  Fig.  21b  by  a  photoluminescence 
current  due  to  beam  excitation  of  the  non-radiating  flow. 

This  behavior  can  be  correlated  with  the  stagnation  pressure 
measurements  in  the  freestream  as  shown  in  Fig.  23,  where 
the  upper  trace  presents  the  "on  beam"  photocurrent  (density) 
measurement  and  the  lower  traces  correspond  to  stagnation 
pressux  e. 

Because  of  the  timewise  difference  in  the  location  of 
the  test  as  obtained  with  the  cylindrical  electrostatic  probes 
and  the  electron  beam  density  probe,  and  the  fact  that  no 
ionized  particles  are  present  in  the  flow  where  density 
measurements  were  made,  it  was  suggested  that  the  density 
measurements  were  actually  made  in  the  region  where  the  driver 
gas.  Helium  in  this  case,  arrives  at  the  test  section.  In 
order  to  verify  that  the  density  measurements  were  made  in 
the  region  where  the  flow  in  the  test  section  consists  of 
the  driven  gas,  air  or  nitrogen  in  our  case,  and  not  Helium, 
the  following  measurements  were  made. 

The  electron  beam  density  probe  was  operated  during 
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tests  conducted  in  the  hypersonic  shock  tunnel  and  the 

photoluminescence  current  recorded  during  each  test.  In  these 

tests  two  different  narrow  band-pass  filters  were  used  in 

the  electro-optical  system  in  order  to  observe  the  0-0  and 

0-1  vibrational  bands  of  the  first  negative  system. 

o  o 

The  one  for  the  0-0  band  was  3914  A  with  100  A  half-band  width 

o  o 

and  the  other  filter  for  the  0-1  band  was  4227  A  with  39  A 

half-band  width.  Since  beam  excited  Helium  radiates  at 

o  o 

3965  A, a  wave  length  which  is  transmitted  by  the  3914/100  A 

filter,  hence  just  using  this  filter  in  these  tests  could  not 

determine  the  arrival  of  the  Helium  at  the  test  section.  In 

order  to  overcome  this  problem,  the  second  interference 
o 

filter  4227/39  A  was  introduced  into  the  electro-optical 

system.  Since  in  this  range  of  wave  length  no  beam  excited 

radiation  of  the  Helium  exists,  hence  the  photoluminescence 

current  is  only  due  to  radiation  of  the  beam  excited  nitrogen. 

By  comparing  the  photocurrent  obtained  during  a  freegtream 

o 

run  with  the  3914/100  A  interference  filter  to  the  photocurrent 

o 

obtained  with  the  4227/39  A  interference  filter,  and  taking 
into  account  the  relative  emission  of  the  0-0  and  0-1 
vibrational  bands,  and  the  transmission  of  each  filter,  it  was 
established  that  the  flow  consists  of  nitrogen. 

B .  Langmuir  Probe  Response  in  Transition  Regime 

An  experimental  investigation  on  the  effect  of  neutral 
density  increase  (mainly  decrease  in  mean  free  path  Xn_n)  on 
the  current  collection  of  cylindrical  probes  in  slightly 
ionized  hypersonic  flow  was  carried  out  during  the  course  of 
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this  work.  This  information  was  needed  in  order  to  reduce 

tho  data  obtained  by  cylindrical  electrostatic  probes  in 

the  flow  field  about  the  plate.  Since  there  exists  a 

compression  zone  at  the  outer  edge  of  the  disturbed  region 

about  the  sharp  leading  edge  flat  plate  in  question  where 

the  local  neutral  density  p  might  be  a  few  times  higher 

than  the  free  stream  density  pm,  it  was  felt  that  the  probes 

used  Rp=0.010"  will  operate  in  a  transition  regime  in  this 

region  and  not  as  free  molecular  cylindrical  probes.  In  order 

to  investigate  the  effect  of  density  increase,  it  was 

decided  to  utilize  the  flow  in  the  nozzle  of  the  HST  as  the 

test  medium.  Here  the  density  could  be  continuously  varied 

maintaining  the  velocity  almost  constant.  The  slight  changes 

in  ion  temperature  Ti  as  one  moves  from  the  test  section  into 

the  nozzle  up  to  about  40  inches  from  the  throat  (length  of 

nozzle  175  inches)  can  be  neglected  particularly  since 
T. 

;~«0  at  any  point  along  the  nozzle  and  the  random  ion 

e 

velocity  as  defined  varies  only  as  the  square  root  of  the 
electron  temperature  Tfi.  Following  the  same  procedure  as  in 
I.ederman's  work  (Ref.  5),  cylindrical  electrostatic  probes  of 
the  same  ?/2Rp  ratio  but  of  variable  diameter  ranging  from 
0.001  inches  to  0.080  inches  were  arranged  in  a  rake  mounted 
on  a  moveable  sting. 

The  different  diameters  of  the  probes  provided  a 
variation  of  R  />  .  in  each  position  where  measurements  were 

p  d 

performed.  All  probes  were  biased  at  -9  volts  (Xp»-35)  for 
ion  collection.  To  correlate  the  collected  ion  current 
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denr.ity  to  the  corresponding  ion  number  density  two  microwave 
resonant  cavities  were  used)  they  were  mounted  on  the  rake 
and  gave  an  independent  measurement  of  the  ion  density  in 
each  test  at  the  different  positions  along  the  nozzle. 

The  sting  with  the  probes  and  cavities  was  moved  stepwise 
from  the  test  section  into  the  nozzle.  In  this  configuration 
the  probes  were  always  aligned  with  the  flow  direction. 

A  number  of  tests  were  run  under  the  same  initial  test 
conditions,  at  different  stations  in  the  nozzle.  The  ion 
densities  as  obtained  by  the  cavities  and  the  ion  current 
of  each  probe  were  recorded. 

The  charged  particle  number  density  obtained  by  the 

microwave  resonant  cavity  measurements  were  compared  to  the 

nonequilibrium  calculations  of  the  nozzle  flow.  Figure  24 

shows  the  experimental  as  well  as  the  theoretical  results. 

The  scc-tter  of  the  experimental  points  is  within  the  accuracy 

one  can  expect  from  the  resonant  cavity  measurements.  The 

electron  number  densities  obtained  by  the  microwave  cavity 

technique  at  each  axial  station  in  the  nozzle  was  used  to 

compute  the  R  /X  parameter  for  each  probe  and  the  corresponding 
p  u 

a  .  The  results  of  these  tests  are  shown  in  Fig.  25  where 
P 

the  parameter  as  a  function  of  is  plotted  for 

two  different  positions  in  the  nozzle.  The  lower  correspond 
to  a  higher  density  as  well  as  a  higher  charge  particle 
density  at  a  station  further  upstream  in  the  nozzle  (see 
Table  I ) . 

One  can  see  from  this  figure  that  within  the  scatter 
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of  the  data  there  seems  to  be  a  decrease  in 


for  the  same 


°V 

R  /'D  when  the  density  increases.  This  result  agrees  in 
principle  with  the  theoretical  prediction  of  Talbot  and  Chou 
(Ref.  30) .  They  have  shown  that  for  cylindrical  Langmuir 
probes  in  the  transition  regime  the  effect  of  collisions 
on  the  probes  is  to  reduce  the  ion  current  to  the  probe  and 
thus  reduce  the  inferred  number  density. 

C .  Ion  Density  Distribution  About  the  Plate 

Cylindrical  electrostatic  probes  arranged  in  a  rake  were 

used  in  the  course  of  this  work  to  obtain  ion  density 

profiles  in  the  field  of  a  sharp  leading  edge  flat  plate 

in  slightly  ionized  rarefied  hypersonic  flow.  Tests  were 

performed  on  the  plate  at  Mach  18  and  zero  angle  of  attack; 

all  tests  reported  in  this  section  were  conducted  in  air 

(commercial  dry  bottled  air),  and  the  flow  conditions  in  the 

test  section  are  listed  in  Table  II.  The  arrangement  of  the 

probe  rake  with  respect  to  the  plate  is  shown  in  Fig.  16. 

As  can  be  seen  from  this  figure,  the  rake  was  perpendicular 

to  the  plate  surface  and  fixed  in  position  with  respect  to 

the  tunnel.  Since  the  rake  was  stationary,  the  plate  was 

moved  in  order  to  obtain  ion  density  profiles  at  various  axial 

positions  along  the  centerline  of  the  plate.  The  probes  on 

the  rake  were  aligned  parallel  to  the  freestream  flow  direction 

and  biased  at  -9  volts  with  respect  to  ground(x  e-35)  for 

P 

ion  collection. 

Figure  26  represents  a  sample  of  the  oscilloscope  traces 
as  obtained  in  these  tests.  As  can  be  seen,  these  traces 
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show,  after  an  initial  peak,  a  region  of  relatively  uniform 
amplitude  which  could  be  considered  as  the  response  in  the 
established  uniform  flow.  This  region  will  be  referred  to 
as  the  flat  portion.  In  reducing  the  data,  only  the  flat 
portion  of  the  probe  response  was  used,  the  reason  being 
that  recently  it  was  shown  by  Bornstein  (Ref.  31  )  that  the 
peak  portion  of  the  probe  response  depends  on  the  load 
resistor  of  the  probe  (R^)  and  the  plasma  properties  in 
question . 

Using  the  procedure  devised  by  Lederman  et  al, (Ref.  5) 
and  discussed  previously,  ion  number  density  was  determined 
from  the  ion  current  collected  by  the  probes.  Since  the  ion 
number  density  depends  on  the  nondimensiortal  current  density 
Op  which  is  a  function  of  itself,  an  iteration  scheme 
had  to  be  used  to  reduce  the  data,  and  a  computer  program 
has  been  devised  which  carried  out  the  data  reduction. 

Figure  27  through  Fig.  30  show  the  ion  density  profiles 
obtained  at  x=3.0,  6.0,  9.0,  12.0,  15.0  inches  along  the  plate. 
The  data  presented  is  in  a  normalized  form,  where  the  ion 
density  *t  each  point  in  a  given  profile  was  calculated  using 
the  flat  portion  of  the  probe  response,  and  then  normalized 
with  the  freestream  value  obtained  in  each  test  by  the  probes 
located  outside  the  disturbed  region  about  the  plate.  In  each 
of  these  figures  is  included  the  calculated  ion  density 
profile  as  obtained  from  the  numerical  solution  of  the  slightly 
ionized  hypersonic  flow  about  the  plate.  In  this  calculation, 
the  following  conditions  and  properties  were  used.  The  free 
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Te 

st  roan,  olcctron  to  ion  temperature  ratio  was  (—•)m  =  42.0.  The 

1  i 

ion  diffusion  coefficient  D.  was  evaluated  as  outlined  in 

1 

Appendix  A.  The  value  of  the  charged  species  concentration 

at  the  wall  used  as  Cw=0,  to  agree  with  the  ion  density 

obtained  near  the  wall  using  the  cylindrical  electrostatic 

probes.  As  was  pointed  out  before,  the  value  of  the  Schmidt 

number  has  a  very  small  effect  on  the  ion  density  profiles 

and  for  the  problem  in  question,  S  =1.0  was  used  in  all  the 

ci 

calculations  carried  out.  The  calculated  ion  density  shows 
very  good  agreement  with  the  measured  values  up  to  the  point 
where  the  outer  compression  region  starts.  At  this  point, 
one  observes  that  the  calculated  ion  density  starts  to  increase 
and  reaches  values  of  one  and  a  half  to  over  three  times  its 
free  stream  value,  something  the  measured  ion  density  profiles 
do  not  exhibit. 

This  may  be  attributed  to  the  fact  that  the  theory  does 
not  apply  at  this  region. 

D .  Density  Profiles  on  the  Plate 

Density  profiles  in  the  disturbed  region  cfoout  a  sharp 
leading  edge  flat  plate  in  a  rarefied  hypersonic  flow  were 
obtained  in  the  course  of  this  work  using  an  electron  beam 
density  probe.  Experiments  were  conducted  in  the  hypersonic 
shock  tunnel  on  a  plate  at  zero  angle  of  attack,  Mach  number 
18,  temperature  ratio  Tw/To=0.075,  and  Reynolds  number  of 
1 'SOO  inch. 

The  following  arrangement  of  the  plate  and  the  optics  was 
used  in  the  density  survey.  The  plate  was  mounted  on  the 
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drift  tube  of  the  electron  gun;  the  beam  was  allowed  to 
emerge  through  a  small  orifice  in  the  plate  surface  and 
impinge  on  the  beam  collector  mounted  12  inches  away  from 
the  plate.  The  radiation  emitted  due  to  the  beam  emerging 
from  the  model  surface  and  traversing  the  test  section  was 
monitored  during  each  test  with  the  electro-optical  system 
which  could  be  positioned  at  any  point  along  the  beam.  The 
surf  ‘ce  position  of  the  optics  was  determined  visually  by 
positioning  the  electro-optical  system  such  that  the  center 
of  the  beam  lens  coincides  with  the  plate  surface. 

Prior  to  each  test  the  electro-optical  system  was 
positioned  at  a  desired  distance  away  from  the  plate  surface 
and  the  photomultiplier  voltage  was  turned  on.  The  electron 
beam  was  in  continuous  operation  before  each  test;  nominal 
operation  conditions  were  20  KeV  accelerating  potential  and 
300  uA  beam  current.  Once  the  test  was  initiated  the  data 
was  displayed  on  Tektronix  565  oscilloscope  and  recorded  by 
a  Polaroid  camera.  In  each  test  the  data  recorded  included 
the  photomultiplier  output  and  the  beam  current  signal.  Since 
ambient  pre-test  pressure  in  the  test  section  was  always 

-4 

less  than  10  mmHg,  and  at  this  pressure  the  photomultiplier 
had  zero  output  (except  for  noise),  hence  the  recorded  portion 
of  the  data  before  the  arrival  of  the  test  gas  gave  the  zero 
signal  level  for  the  recorded  data.  Figure  31  presents  typical 
data  traces  recorded  during  a  test.  The  upper  trace  corresponds 
to  the  photomultiplier  output.  As  can  be  seen,  after  some 
initial  disturbance  there  is  a  gradual  increase  in  the 
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photoluminescence  current  which , after  a  certain  period  of  time, 
levels  off  to  a  constant  value.  The  lower  trace  corresponds 
to  the  beam  current  signal,  which  remains  relatively  constant 
throughout  the  test. 

A  series  of  density  profiles  normal  to  the  plate  surface 
were  obtained,  ranging  from  2  to  9  inches  from  the  leading 
edge.  Because  of  the  extremely  short  test  duration  pointwise 
measurements  had  to  be  performed  in  order  to  obtain  a  density 
profile.  Figure  32  through  Fig.  35  present  the  normalized 
density  profiles  in  the  disturbed  region  about  the  plate.  The 
density  at  each  point  of  the  profile  was  normalized  with 
respect  to  the  freestream  density  after  corrections  were 
made  for  changes  in  beam  current  between  tests  and  taking 
into  account  variation  in  the  photomultiplier  output  depending 
on  its  position  with  respect  to  the  plate.  It  should  be 
noted  here  that  the  accuracy  of  positioning  the  electro-optical 
system  at  any  given  point  normal  to  the  plate  was  ±0.025  inches. 
This  fact  had  an  effect  on  the  accuracy  of  the  density 
measurements  in  the  shock  layer  where  rapid  density  changes 
occur  in  a  relatively  narrow  region. 

E*  Flush  Probe  Measurements  and  Size  Study 

Flush  probe  ion  saturation  current  collection  gives 
information  concerning  the  ambipolar  diffusion  flux  to  the 
plate  surface  and  the  charge  number  density  at  the  sheath 
edge  of  the  plate.  Via  boundary  layer  theory,  charged 
particle  density  at  the  boundary  layer  edge  can  also  be 
determined.  In  this  experimental  investigation,  flush 
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probes  were  used  only  to  determine  the  ambipolar  diffusion 
flux  to  the  plate  surface. 

Measurements  were  also  made  of  the  current  voltage 
characteristics  of  the  flush  probes.  This  was  achieved  by 
changing  the  bias  of  the  different  diameter  probes  at  each 
test  a;.d  recording  the  collected  ion  current.  This  pointwise 
measurement  of  the  current  voltage  characteristic  was  due 
to  the  fact  of  extremely  short  test  duration  whereby  no 
sweeping  of  the  probe  bias  was  possible  during  a  test. 

Figure  39  shows  a  typical  current  voltage  characteristic  of 
a  circular  flush  probe  biased  as  a  single  probe.  For  all 
probes  large  and  small,  except  the  3/4"  diameter,  no  ion 
saturation  current  could  be  observed.  Even  for  the  largest 
probe,  3/4”,  the  ion  saturation  was  ambiguous.  The  ion 
current  was  found  to  be  a  very  strong  function  of  the  bias 
potential  and  probe  size.  Therefore,  a  series  of  size  studies 
was  undertaken  to  investigate  these  effects. 

Flush  probe  studies  were  made  by  using  the  ion  current 

collected  by  the  3/4"  diameter  probe  at  a  bias  of  >3  volts, 

as  a  reference.  The  reason  for  this  choice  was  the  fact  that 

since  no  ion  saturation  current  was  reached  even  for  this 

probe  as  shewn  in  Fig.  40,  the  low  value  of  would  partially 

satisfy  the  desire  to  minimize  the  sheath  size  which  is 

3/4 

proportional  to  V  and  at  the  same  time,  provide  a 
reasonable  current  signal.  Thus,  by  taking  the  ion  current 
collected  by  the  3/4"  flush  probe  at  a  bias  of  -3  volts, 
it  was  felt  that  fringing  effects  were  kept  minimal  and  that 
this  current  is  close  to  the  ion  saturation  current.  Ion 
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current  density  Ji#  to  a  smaller  sized  flush  probe  at  a 
particular  x  was  then  normalized  by  the  current  density  of 
a  3/4"  probe.  This  normalized  current  density  was  then 
plotted  against  the  parameter  dg/Rp  as  shown  in  Pig.  41  to 
demonstrate  the  fringing  field  effects;  here  dg  is  the 
sheath  thickness  at  the  biased  flush  probe  and  the  probe 

radius.  For  the  determination  of  dg  use  was  made  here 
of  the  simple  free-fall  sheath  relationship  given  by  the 
Child-Langmuir  Law  cited  by  Brown  (Ref.  32), 

*  ft  df*  (23) 

where 

<*  -% 

XD  is  the  Debye  shielding  distance  based  on  the  plasma 
conditions  at  the  sheath  edge.  Xp  was  calculated  from  the 
ion  density  near  the  surface  of  the  plate  as  measured  by 
the  cylindrical  electrostatic  probes, and  the  electron 
temperature  Te  was  assumed  to  be  constant  and  equal  to  its 
free st ream  value  since  the  flow  is  frozen. 

The  combined  sheath  area  effect  was  first  demonstrated 
by  applying  a  constant  bias  potential  to  all  probes  of 
different  sizes  under  the  same  flow  conditions.  To  vary 
d  /R  even  more,  different  bias  potentials  up  to  -200  volts 

a 

were  applied  to  these  probes. 

F.  Ambipolar  Diffusion  Flux  Along  the  Flat  Plate 

The  ambipolar  diffusion  flux  of  the  ionized  particles 
to  the  plate  surface  at  different  axial  locations  along  the 
plate  was  obtained  during  the  course  of  this  work  using  three 
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different  approaches. 

In  the  first  approach,  following  the  formation  of  Chung 
(Ref.  16)  the  current  density  at  the  plate  surface  is  related 
to  the  ambipolar  diffusion  coefficient  and  the  gradient  of 
ion  number  density  at  the  plate  surface  in  the  following 
manner  s 

(fsMw  (24) 

where 

°a  *  ^  (l  ♦  ^-) 

Since  of  interest  was  the  ambipolar  diffusion  flux  along  the 
plate,  hence,  the  ion  current  density  distribution  as  a 
function  of  x  was  calculated  using  this  relation  (Eq.  24) 
together  with  the  numerical  calculation  of  the  slightly  ionized 
hypersonic  flow  about  the  plate  in  question.  It  was  assumed 
here  that  the  ion  and  neutral  particle  temperature  are  equal 
(T^-T)  and  that  the  electron  temperature  is  constant 
throughout  the  region  since  the  flow  is  frozen  (Te»3000°K) . 

The  ambipolar  diffusion  coefficient  was  determined  at  each 
point  as  outlined  in  Appendix  A. 

The  second  method  consisted  again  of  utilizing  Eq.  (24) 
together  with  the  value  of  at^/dyfw  obtained  from  the  experi¬ 
mental  ion  density  profiles  at  the  various  axial  station. 

The  density  and  temperature  distribution  in  the  disturbed 
region  about  the  plate  was  obtained  here  from  the  numerical 
calculation . 

The  third  method  used  to  establish  the  ambipolar  diffusion 
flux  to  the  plate  was  to  measure  the  ion  current  collected  by 
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the  3/4"  diameter  flush  electrostatic  probes  imbedded  at 
different  axial  positions  along  the  plate  and  biased  at  -3 
volts  for  ion  collection. 

Figure  42  presents  the  ambipolar  diffusion  flux  along 
the  surface  of  the  plate  as  obtained  using  those  three  methods. 
Figure  43  shows  the  non-dimensional  current  density 
distribution  along  the  plate  where  the  current  density 
collected  by  the  flush  probes  was  non-dimensionalized  with 
the  current  collected  by  the  probe  located  at  the  furthest 
downstream  position  (x=16h  inches)  from  the  leading  edge. 
Included  in  this  figure  is  the  non-dimensional  ambipolar 
flux  of  the  ionized  particles  as  obtained  from  the  numerical 
solution  for  the  flow  in  question. 

G.  Cone  Wake  Density  Measurements 

As  mentioned  above  in  previous  work  by  Lederman  and 
Avidor  (Ref.  7), it  was  suggested  that  ions  in  a  preionized 
flow  can  be  used  for  flow  field  tracing  and  that  in  chemically 
frozen  flow  fields,  the  ion  density  distribution  measured  can 
be  interpreted  in  terms  of  neutral  density  distribution, 
provided  diffusion  is  negligible.  In  this  work  and  another 
investigation  (Ref.  8) , ion  density  distribution  in  the  near 
wake  of  slender  and  blunted  bodies  was  obtained  in  a  slightly 
preionized  hypersonic  flow  generated  in  a  hypersonic  shock 
tunnel.  The  present  investigation  of  density  measurements 
in  the  near  wake  of  a  cone  was  initiated  because  the  wake 
provides  the  different  flow  fields  in  which  diffusion  is  more 
or  less  significant,  and  thus  provides  the  necessary  flow  field 
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to  determine  the  correctness  or  incorrectness  of  the  flow  field 
diagnostics  as  suggested  by  Lodorman  et  al. 

Utilizing  an  electron  beam  density  probe  for  this 
purpose,  tests  were  conducted  on  a  sharp  10°  half  angle 
cone  with  a  4-inch  base  diameter  at  Mach  number  18  and  zero 
angle  of  attack.  All  tests  on  the  cone  wake  reported  here 
were  conducted  in  pure  nitrogen  (commercial  dry  bottled 
nitrogen),  and  the  flow  conditions  in  the  tost  section  are 
listed  in  Table  11. 

For  any  wake  measurements,  the  model  suspension  system 
is  of  critical  importance.  For  the  present  tests,  the  model 
was  suspended  on  0.020  inch  diameter  wire.  It  was  found  in 
previous  work  by  Lederman  (Ref.  29)  that  this  size  wire  did 
not  introduce  any  observable  disturbance  in  the  wake  region 
under  the  same  flow  conditions  as  those  in  question  here. 

A  further  experimental  detail  concerned  the  positioning 
of  the  model  with  respect  to  the  electron  beam.  The  model 
was  suspended  on  wires  from  a  frame  mounted  to  a  movable 
sting,  which  allowed  one  to  position  the  cone  at  different 
axial  stations  (x/Dn)  with  respect  to  the  beam's  position  which 
was  fixed.  The  x/D  =0  position  was  obtained  before  each  test 
by  projecting  the  beam  through  the  test  section  and  letting 
it  impinge  on  the  shoulder  of  the  cone,  and  from  this 
reference  point,  the  cone  was  moved  to  a  desired  x/Dn  position 
by  moving  the  sting.  Thus,  after  the  cone  was  set  in  position 
at  a  given  x/D^#  the  electro-optical  system  was  moved  to  a 
desired  radial  position  by  moving  it  along  the  beam.  In  each 
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test  two  data  points  were  obtained  in  any  radial  profiles#  each 
by  one  electro-optical  system.  Because  of  the  short  test 
duration,  one  cannot  move  the  electro-optical  systems  to 
obtain  a  continuous  density  profile  during  the  test  run. 

Density  measurements  obtained  in  the  near  wake  of  the  cone 
include  axial  density  distribution  from  x/DB»0.25  to  x/bB»  3.0  . 
At  x/Db*1.0  and  x/DB«3.0  a  radial  density  profile  was  obtained 
which  extended  to  the  outer  edges  of  the  flow  field.  Before 
presenting  the  data,  one  aspect  of  the  work  should  be 
emphasized.  All  of  the  density  measurements  are  direct 
"point"  measurements.  The  resolution  of  "point"  size  in 
this  case  is  2  mm.  Figure  44  shows  an  oscilloscope  trace 
of  the  electroluminescence  intensity  for  a  freestream  run 
and  one  taken  in  the  wake  of  the  cone.  The  traces  include 
a  sharp  change  in  photocurrenc  occurring  when  the  flow 
arrives  at  the  beam  station  and  then  the  photocurrent  increases 
for  awhile  and  levels  off  to  a  constant  value  over  a  period 
of  0. 5-1.0  msec.  The  traces  also  include  the  beam  current 
as  monitored  during  each  test. 

The  fluctuations  observed  in  the  current  trace  are 
due  to  the  fact  that  the  beam  collector  which  was  biased  at 
+9  volts  with  respect  to  ground  was  actually  a  biased 
electrostatic  probe  collecting  electrons  from  the  slightly 
ionized  flow.  This  current,  however,  was  found  to  be  negligible 
in  the  region  where  the  tests  were  conducted.  It  was  found 
that  no  observable  attenuation  of  the  beam  current  did  exist 
during  the  test  duration.  It  is  also  worth  mentioning  here 
that  during  the  course  of  this  work,  it  was  found  that  the 


repeatability  of  tests  was  excellent. 

Figure  45  presents  the  axial  density  distribution  in 
the  near  wake  of  the  cone  as  obtained  in  the  present  work. 

The  density  at  each  point  on  the  axis  of  the  wake  was 
normalized  with  the  freestream  density.  Included  in  this 
figure  for  the  purpose  of  comparison  is  the  ion  density 
distribution  in  the  near  wake  of  the  cone  as  obtained  by 
Lederman  et  al.*with  cylindrical  electrostatic  probes  (Ref.  29). 
The  next  two  figures.  Fig.  46  and  47.  show  two  radial 
profiles  measured  at  x/DB*1.0  and  x/Dga3.0.  Those  were 
obtained  by  positioning  the  cone  at  the  particular  x/Dfi  and 
measuring  the  density  by  moving  the  electro-optical 
system  in  the  radial  direction,  starting  at  the  wake  axis 
and  mewing  out  to  the  outer  edges  of  the  flow  field.  Again, 
the  density  is  presented  in  a  normalized  form  where  at  each 
point  the  density  was  normalized  with  the  local  freestream 
value.  These  figures  also  include  the  radial  ion  density 
profiles  as  obtained  by  Lederman  et  al. (Ref .  29). 
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VII.  DISCUSSION  OF  RESULTS 


As  pointed  out  above  and  shown  in  Fig.  24  through 
Fig.  27,  the  normalized  profiles  of  the  ion  density  about  the 
flat  plate  at  the  several  stations  investigated  agree  within 
the  experimental  error  remarkably  well  with  the  theoretically 
calculated  profiles  up  to  the  point  where  the  outer 
compression  region  starts.  At  this  point  one  can  observe  a 
departure  between  the  measured  and  calculated  profiles. 

Although  there  is  no  complete  explanation  of  this  behavior 
possible  at  this  time,  it  should  be  noted  that  the  absolute 
ion  density  values  may  have  been  underestimated  through  the 
use  of  the  upper  curve  in  Fig.  25.  This  curve  was  obtained 
at  a  density  which  was  lower  than  the  actual  density  in  the 
compression  region,  as  seen  in  the  same  figure,  the  lower  curve- 
an  increase  in  density  would  result  in  a  decrease  in 
which  would  in  turn  result  in  a  higher  number  density  value 
at  the  given  point.  Furthermore,  it  should  be  indicated  that 
the  numerical  solution  in  the  compression  region  is  of 
questionable  validity. 

In  Figs.  32  through  35,  the  experimentally  obtained 
neutral  density  profiles  normal  to  the  plate  are  shown. 
Superimposed  on  these  are  the  theoretically  obtained  profiles. 

It  is  evident  from  these  figures  that  good  agreement  between 
the  theoretical  and  experimental  density  profiles  exist. 

Slight  deviations  can  be  observed  in  the  position  of  the 
maximum  density.  The  absolute  values  of  the  theoretical 
and  experimental  maximum  density  (pfoax/Pe„)  deviate  more 


50 


significantly.  This  deviation  appears  to  decrease  with  the 

increase  in  x  (decrease  in  Vw  ) (see  Fig.  36).  This  kind  of 

x 

behavior  was  observed  by  a  number  of  other  investigators 

(Refs.  28/33,34).  The  next  two  figures.  Figs.  37  and  38. 

consist  of  a  superposition  of  an  experimentally  obtained 

neutral  and  ion  density  profile.  Because  of  the  large 

difference  in  Reynolds  number  for  the  ionized  flow  and  the 

nonionized  flow  (see  Table  II)  the  correlating  parameter  used 

in  this  superposition  was  Vn  . 

x 

As  can  be  seen,  there  is  no  agreement  between  the  two 

profiles.  In  the  viscous  region  the  ion  density  drops  down 

to  much  lower  values  than  the  density.  This  could  partly 

be  explained  by  the  large  ambipolar  diffusion  of  the  ionized 

particles  because  of  the  relatively  large  electron  to  ion 

Te 

temperature  ratio  (=—)  under  the  test  conditions  in  question. 

Ti 

In  the  compression  region  the  density  shows  a  definite 
overshoot  whereas  the  ion  density  does  not. 

In  Fig.  42, the  ambipolar  diffusion  flux  of  the  ionized 
particles  to  the  plate, as  determined  experimentally  and 
theoretically,  is  shown.  It  is  clearly  evident  from  this 
figure  that  the  theoretical  and  experimental  values  may 
deviate  by  a  factor  of  two.  This  deviation  could  be  ascribed 
partially  to  the  experimental  scatter  of  the  data,  the 
fringing  effects  on  ion  current  collection  by  the  flush 
electrostatic  probes  as  indicated  in  Fig.  41,  and  the 
uncertainty  in  the  calculated  values  because  of  initial 
condition  effects  on  the  numerical  solution  for  Vm  >  0  .4 

cxj 
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( x  «r  4)  .  Figure  43  presents  the  nondimensionlized  ambipolar 
flux  to  the  plate.  Again,  significant  deviation  occurs  for 
x<  4  . 


In  the  next  few  figures  the  comparison  between  ion 
and  neutral  density  distributions  in  the  near  wake  of  the 
cone  is  shown.  Figure  45  presents  the  axial  distribution  and 
as  one  can  see,  there  is  a  significant  deviation  between 
the  ion  and  neutral  density  up  to  about  x/Dfi=2.0.  Figure  46 
shows  the  radial  distribution  and  as  is  evident  from  this 
figure,  the  difference  between  the  ion  and  neutral  density 
is  in  particularly  large  in  the  core  of  the  wake.  It  should 
be  noted  that  the  overshoot  in  the  density  and  ion  density 


indicating  the  location  of  the  shock  is  at  the  same  radial 

• .  •  .  _ . _ .  .  i  _  i _ .  *  /  ^nax  i _ » 


position  and  about  the  same  strength  ( 


max  ) . 


Figure  47  shows  the  radial  distribution  at  x/b  «3.0j  in  this 


case  very  good  agreement  exists  between  the  ion  and  neutral 


density  profiles  even  in  the  center  core  region. 


In  view  of  the  above,  it  is  quite  evident  that  the 


electrostatic  probe  technique  can  only  be  utilized  for  flow  field 


determination  in  regions  where  diffusion  of  the  ionized 
particles  is  negligible.  In  these  regions  a  one-to-one 
correlation  can  be  obtained  for  the  ion  and  neutral  densities. 


Finally,  it  should  be  noted  that  two  flow  regimes  were 
utilized  in  the  course  of  this  investigation.  One,  containing 
ionized  particles  of  approximately  0.4  msec,  duration;  the 
other  of  a  higher  density  following  the  ionized  flow  region 
and  containing  a  negligible  amount  of  ionized  particles. 
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VIII.  SUMMARY  AND  CONCLUSIONS 


Two  diagnostic  techniques,  the  cylindrical  electrostatic 
probe  and  the  electron  beam  density  probe,  were  applied  to 
flow  field  surveys  in  the  disturbed  region  aoout  a  flat 
plate  and  in  the  near  wake  of  a  10°  half  angle  cone  in  a 
rarefied  slightly  ionized  hypersonic  flow.  Their  relative 
merits  have  been  evaluated.  It  was  found  that  both  can  be 
applied  for  flow  field  investigations  and  particularly  density 
tracing  and  mapping,  provided  that  diffusion  effects  are 
nonexistent  or  negligible.  In  canes  where  diffusion  is 
significant,  the  electrostatic  probe  technique  may  provide 
highly  misleading  data.  The  correspondence  between  the  ion 
and  neutral  densities  in  this  case  are, on  a  one-to-one  basis, 
nonexistent.  On  the  other  hand,  the  electron  beam  density 
probe  is  applicable  in  diffusive  and  nondiffusive  regions. 

Its  limitations,  however,  are  at  extremely  low  densities 
where  the  signal -to-noise  ratio  becomes  intolerable,  or  at 
the  high  density  range  where  collisions  cause  quenching  of 
the  beam.  From  the  above,  it  is  evident  that  both  of  those 
techniques  are  supplementary}  whereas  the  former  has  no 
limitation  on  the  density  regime,  the  latter  is  not  limited 
by  diffusion  effects. 

It  can,  therefore,  be  concluded  that: 

i)  The  electrostatic  probe  can  be  utilized  as  a  field 
tracing  tool  in  regions  where  diffusion  is  nonexistent  or 
negligible  irrespective  of  density. 

ii)  The  electron  beam  density  probe  provides  a  unique 
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technique  for  tracing  and  determination  of  flow 
configurations  irrespective  of  diffusion  effects,  provided 
it  is  operated  in  the  proper  density  regime. 

iii)  Ion  density  profiles  normal  to  the  plate  agree  very 
well  in  the  viscous  layer  with  results  obtained  from  the 
numerical  solution  using  the  Rudman-Rubin  theory  for  the 
hypersonic  flow  about  the  plate  modified  by  Chung's  formulation 
for  the  ambipolar  diffusion  of  the  ionized  particles  in  the 
disturbed  region. 

iv)  Neutral  density  profiles  obtained  normal  to  the  plate 
indicate  a  certain  deviation  from  the  theoretically  predicted 
result.  The  experimentally  obtained  data  agree  with  data 
obtained  by  other  investigators  under  similar  flow  conditions. 
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APPENDIX  A 


DIFFUSION  AND  AMBI POLAR  DIFFUSION  COEFFICIENT  FOR  No+  IN  AIR 

Most  data  on  the  motion  of  ions  in  neutral  gases  is 
given  in  terms  of  mobility,  which  is  related  to  the 
diffusion  coefficient,  D^#  by  Einstein's  relation 

Pl  /i  (A.l) 

the  mobility  which  is  usually  given  in  literature 
represents  the  value  of  at  the  standard  gas  number  density 
of  Nq«2.69  x  10® 8  particles/m3,  which  corresponds  to  standard 
conditions  of  760  Torr  and  273°k.  The  mobility  at  other 
densities  is  related  to  by  the  equation 

/* l  (N,rj-j£-/i.(T)  <A.2> 

but  NQ/to  *  pQ/p,  pQ  is  the  gaa  density  at  standard  conditions. 

Mason  (Ref.  35)  gives  values  for  uiQ(T)  for  N0*  in 
air  which  can  be  approximate  d  by  the  following  expression 

(rj«  2.73  «■  1.466  «  10’*T  [■Sslttr]  (A. 3) 

in  the  range  0°k  <  T  <  1000°k.  The  asfeipolar  diffusion 
coefficient  is  related  to  the  diffusion  coefficient  of  the 
ion  in  the  following  ways 

0L  (lef-)  (A. 4) 

Substituting  Eq.  (A.l)  into  (A. 4)  and  using  the  expression 
(N,T) ,  we  obtain 

D*  -  ~  (t(  «■*)  75  ♦1.444*  icT*Ti]  (A. 5) 
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Eq.  (A. 5)  was  used  in  all  the  calculations  concerning  ion 
mobility  in  the  flow. 
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TABLE  I 


FLOW  CONDITIONS  ALONG  THE  NOZZLE 
(IONIZED  FLOW  REGION) 


X 

40  in. 

Hi 

1 

1 

T(°k) 

315 

225 

167 

115 

P(Atm) 

3.96  x  10"4 

SBSSEB 

10 

p(Atm) 

4.75  x  10"4 

[  5 

9.4  x  10" 

EBHRI 1 

,  mom 

1.5  x  1010 

9  x  109 

4  x  10® 

Mi 

Te(°k) 

3000 

3000 

3000 

msi _ 
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TABLE  II 


FLOW  CONDITIONS  IN  THE  TEST  SECTION 

Ionized  Flow  Region  Nonionized  Flew  Region 


T 

00 

~70°k 

6 

~70°k 

—  6 

2.6  x  10"  Atm 

i*3  x  10  Atm 

^cr> 

1.1  x  10-5Atm 

5. 5  x  10  Atm 

M 

00 

~18.0 

~18 . 0 

R 

eoo 

~300/in . 

~1500/in . 

Ti„ 

-70°k 

T 

e 

00 

~3000°k 

N. 

1co 

3  x  10a  el/cm3 

Vn 

~0.1  in. 

~0.020  in. 

0.01  in. 
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FIG.  I  FLOW  PATTERN  ABOUT  SHARP  LEADING  EDGE  FLAT  PLATE 
M  HYPERSONIC  FLOW 


SOLUTION  VALIDITY 
QUESTIONABLE  IN  THIS 
REGION  y 


0.4  0.8 


1.6  2.0 


W 


FIG.  2  EFFECT  OF  Cw  ON  CALCULATED  ION 
DENSITY  PROFILE  ,, 


0  0.4  0.8  1.2  1.8  2.0  24  2.8 

FIG.  3  EFFECT  OF  SCHMIDT  NUMBER  ON 

CALCULATED  ION  DENSITY  PROFILE 
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FIG.  4  KM  CURRENT  VS  X0/R,  FOR  VARIOUS 

VALUES  OF  -ty/kV  CYLINDRICAL  PROBE; 
T|  /\  -0. 


10*  2  4  •  •  to"'  2  4  6  *  10° 


FIG.  5  DIMENSIONLESS  ION  CURRENT  aP  AS 
FUNCTION  OF  Rp/X0 


4.0 


Xp-35 


//2Rp 

FIG.  6  NONDIMENSION  AL  CURRENT  DENSITY  AS  A  FUNGTIOI* 
OF  THE  PROBE  LENGTH  TO  DIAMETER  RATIO 
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FIG.  8 
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nc  9  «OOUTIC  DIAGRAM  OF  ELECTRON  BEAM 


ELECTRON  GUN  CONTROL  SYSTEM  ’ 


X  (A) 


X  (A) 


FIG. 13  THE  EMI  PHOTOMULTIPLIER  QUANTUM 

EFFICIENCY  AND  BELOW  THE  TRANSMISSION 
CURVE  OF  THE  INTERFERENCE  FILTER 


* 


2 
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FIG.  14  THE  ELECTRO- OPTICAL  SYSTEM  AND  PHOTOMULTIPLIER  CIRCUITRY 
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R6. 15  LANGMUIR  PROBE  CONSTRUCTION  AND 
CIRCUIT 
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FIG.  17  SETUP  FDR  DENSITY  MEASUREMENTS  ABOUT  THE  PLATE 


FIG.  18  PHOTOGRAPHIC  VIEW  OF  THE  PLATE, 
ELECTRON  GUN  AND  ELECTRO -OPTICAL 
SYSTEMS 
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FIG.  19  EXPERIMENTAL  SETUP  FOR  CONE  WAKE  DENSITY  SURVEY 
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FI6.2I  (a)  RESPONSE  OF  CYLINDRICAL 
ELECTROSTATIC  PROSES 
(b)  FREE  STREAM  PHOTOCURRENT, 
NO  ELECTRON  SEAM 
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FIG.  22  (a)  FREE  STREAM  PHOTOCURRENT,  NO 
ELECTRON  BEAM 

(b)  FREE  STREAM  PHOTOCURRENT,  WITH 
ELECTRON  BEAM 


PHOTOCURRENT 


(o) 

BEAM  CURRENT 


(b) 

STAGNATION 

PRESSURE 


hlH— 2  - H 

1.  IONIZED  FLOW  REGION 

2.  NONIONIZED  FLOW  REGION 


FIG.  23  (o)  FREE  STREAM  PHOTOGURRENT,  WITH 
BEAM 

(b)  FREE  STREAM  STAGNATION  PRESSURE 
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ELECTRON  DENSITY  (el/cm3) 


FIG.24  ELECTRON  NUMBER  DENSITY  ALONG  THE  NOZZLE 
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FIG.  25  DIMENSIONLESS  ION  CURRENT  (ap)  DEPENDENC 
ON  NEUTRAL  MEAN  PATH  Xn.n 
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FIG.  26  RESPONSE  OF  ELECTROSTATIC  PROBES  IN 
THE  FLOW  ABOUT  THE  PLATE.  TOP  TO 
BOTTOM  REPRESENT  PLATE  SURFACE  TO 
FREE  STREAM  RESPONSE 
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X  *3.0"  X  *6i0 

•  EXPERIMENT  •  EXPERIMI 


NORMALIZED  ION  DENSITY  PROFILE  NORMAL  TO  THE  PLATE 


FIG.  28  NORMALIZED  ION  DENSITY  PROFILE 
NORMAL  TO  THE  PLATE  AT  X-9.0* 
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FI6.2I  (a)  RESPONSE  OF  CYLINDRICAL 
ELECTROSTATIC  PROSES 
(b)  FREE  STREAM  PHOTOCURRENT, 
NO  ELECTRON  SEAM 
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FIG.  22  (a)  FREE  STREAM  PHOTOCURRENT,  NO 
ELECTRON  BEAM 

(b)  FREE  STREAM  PHOTOCURRENT,  WITH 
ELECTRON  BEAM 
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1.  IONIZED  FLOW  REGION 

2.  NONIONIZED  FLOW  REGION 


FIG.  23  (o)  FREE  STREAM  PHOTOGURRENT,  WITH 
BEAM 

(b)  FREE  STREAM  STAGNATION  PRESSURE 
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ELECTRON  DENSITY  (el/cm3) 


FIG.24  ELECTRON  NUMBER  DENSITY  ALONG  THE  NOZZLE 
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FIG.  25  DIMENSIONLESS  ION  CURRENT  (ap)  DEPENDENC 
ON  NEUTRAL  MEAN  PATH  Xn.n 
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FIG.  26  RESPONSE  OF  ELECTROSTATIC  PROBES  IN 
THE  FLOW  ABOUT  THE  PLATE.  TOP  TO 
BOTTOM  REPRESENT  PLATE  SURFACE  TO 
FREE  STREAM  RESPONSE 
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X  *3.0"  X  *6i0 

•  EXPERIMENT  •  EXPERIMI 


NORMALIZED  ION  DENSITY  PROFILE  NORMAL  TO  THE  PLATE 


FIG.  28  NORMALIZED  ION  DENSITY  PROFILE 
NORMAL  TO  THE  PLATE  AT  X-9.0* 


FIG.  29  NORMALIZED  ION  DENSITY  PROFILE 
NORMAL  TO  THE  PLATE  AT  X  ■  12.0" 
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FIG.  30  NORMALIZED  ION  DENSITY  PROFILE  NORMAL 
TO  THE  PLATE  AT  X  - 15.0* 
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FREE  STREAM 


FIG.  31  PHOTOCURRENT  OURlNG  A  TEST 
(BEAM  ON) 
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RG.33  NEUTRAL  DENSITY  NORMAL  TO  THE  PLATE  AT  X»3 


FIG.  34  NEUTRAL  DENSITY  NORMAL  TO  THE  PLATE  AT 
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THEORY  (RUDMAN  ANO 
RUBIN) 


FIG.  35  NEUTRAL  DENSITY  NORMAL  TO  THE  PLATE  AT 
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FIG.  36  DENSITY  RATI 
WAVE  IN  THE 
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FIG.  37  NORMALIZED  ION  AND  NEUTRAL  OENSITY 
PROFILES  NORMAL  TO  THE  PLATE, 

Ve^O.28 


FIG.  38  NORMALIZED  ION  AND  NEUTRAL  DENSITY 
PROFILES  NORMAL  TO  THE  PLATE, 

VflDv»  0.23 
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FIG.  40 


J^-CUARENT  OENSITY  OF  3/4  M.  DIA. 
PROBE  AT  ^>-3v 
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FIG  .  41  NON  DIMENSION  AL  CURRENT  DENSITY  AS  A 
FUNCTION  OF  THE  RATIO  OF  Rp/ds 


RG.  42  AMBIPOLAR  DIFFUSION  FLUX  TO  THE  PLATE 


CURRENT  DENSITY  COLLECTED 


FIG.  43  NONOIMENSIONAL  AMBIPOLAR  FLUX  TO  THE 
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RG.  44  PHOTOCURRENT  TRACES  DURING 
A  TEST  (BEAM  ON) 
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FIG.  45  NORMALIZED  AXIAL  NEUTRAL  DENSITY  AND 
ION  DENSITY  DISTRIBUTION  IN  THE  NEAR 
WAKE  OF  THE  CONE 
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FIG. 46  NORMALIZED  RADIAL  NEUTRAL  AND  ION 
OENSITY  PROFILE  IN  THE  RAKE  OF  THE 
CONE  AT  X/Do  -1.0 
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FIG.  47  NORMALIZED  RADIAL  NEUTRAL  AND  ION 
DENSITY  PROFILE  IN  THE  WAKE  OF  THE 
CONE  AT  X/Dg-3.0 
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